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We discuss a method to create two-dimensional traps as well as atomic shell, or bubble, states 
for a Bose-Einstein condensate initially prepared in a conventional magnetic trap. The scheme 
relies on the use of time-dependent, radio frequency-induced adiabatic potentials. These are shown 
to form a versatile and robust tool to generate novel trapping potentials. Our shell states take the 
form of thin, highly stable matter-wave bubbles and can serve as stepping-stones to prepare atoms in 
highly-excited trap eigenstates or to study 'collapse and revival phenomena'. Their creation requires 
gravitational effects to be compensated by applying additional optical dipole potentials. However, in 
our scheme gravitation can also be exploited to provide a route to two-dimensional atom trapping. 
We demonstrate the loading process for such a trap and examine experimental conditions under 
which a 2D condensate may be prepared. 
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I. INTRODUCTION 

Due to the rapid advances in experimental and theo- 
retical atom optics, it has now become possible to cool 
atoms down to extremely low temperatures. An impor- 
tant characteristic of such ultracold atoms is their sensi- 
tivity to very weak external potentials. In the context of 
atomic Bose-Einstein condensation, this feature is par- 
ticularly attractive. The condensate can be described, to 
a very good degree of approximation, in terms of a sin- 
gle macroscopic wave function which in this way can be 
subjected to intricate probing and manipulation. Conse- 
quently, in recent years various ingenious techniques for 
handling Bose-Einstein condensates (BECs) have been 
developed that exploit this sensitivity, e.g., magnetic 
quadrupole and TOP traps, shallow optical dipole traps, 
phase-imprinting methods to create solitons or vortices, 
and radio-frequency (RF) output coupling, to name but a 
few prominent examples (for a review of these and other 
techniques, see [1]). However, the development of further 
tools and methods still presents a significant objective of 
current work, and forms the basis for advancing research 
in several main areas of BEC physics. 

One such area concerns the experimental realization 
of low-dimensional Bose-Einstein condensates of trapped 
dilute atomic gases. Under reduced dimensionality, the 
condensate properties differ drastically from the well- 
studied three-dimensional case and have been under in- 
tense theoretical debate for some time [2, 3]. Low- 
dimensional BECs are characterized by the fact that due 
to strong confinement by the external trapping poten- 
tial one or more motional degrees of freedom become 
quantum-mechanically frozen before the condensation 
point is reached. Major obstacles, that had been con- 
fronting the preparation of one- or two-dimensional BECs 
for a long time, concerned the construction of suitable 



traps — e.g., very high field gradients would be required 
for magnetic traps — as well as the development of effi- 
cient loading procedures. Only recently has it become 
possible to overcome these difficulties, and the first real- 
izations of atomic BECs in two dimensions have been re- 
ported [4, 5, 6]. In these experiments, the dimensionality 
is reduced by exposing the atoms to steep optical poten- 
tials. In Ref. [4], single one- or two-dimensional BECs 
are prepared, whereas in the experiments of Refs. [5] and 
[6] the use of optical lattices led to the creation of arrays 
of two-dimensional condensates. However, in addition 
to these works, various other interesting ideas regarding 
the manufacture of two-dimensional traps have been put 
forward over the last few years. In these proposals, trap- 
ping is provided either by optical [7] or by magnetic [8] 
means, whereas loading is accomplished through optical 
pumping. 

The purpose of the present paper is twofold. On the 
one hand we wish to promote field-induced adiabatic po- 
tentials as a new versatile tool to manipulate ultracold 
atoms and, in particular, Bose-Einstein condensates. On 
the other hand, we show that they offer a novel route to 
the creation of two-dimensional traps for BECs. An ini- 
tial outline of our results was given in Ref. [9] . Adiabatic 
potentials arise whenever two or more internal atomic 
states, that experience different potentials for the atomic 
center-of-mass motion, are coupled by a strong resonant 
external field. The atomic motion is then no longer dom- 
inated by the different 'bare' potentials but is usefully 
described in terms of the adiabatic potentials, that arise 
from the diagonalization of the bare potentials and the 
couplings at each spatial point. Adiabatic potentials have 
been in use for some time as the underlying mechanism 
for evaporative cooling. However, to our knowledge, lit- 
tle attention has been paid to the fact that they offer a 
lot of further possibilities to control quantum-mechanical 
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atomic motion (one fairly recent application is described 
in [10]). In this paper we will give some examples of the 
application of adiabatic potentials as a tool for manip- 
ulating matter waves, and we hope to stimulate further 
research in this direction. 

After giving a qualitative introduction to adiabatic po- 
tentials in Sec. II, we continue in Sec. Ill to show how 
they can be employed to create so-called 'matter-wave 
bubbles' from a BEC initially trapped in the ground state 
of a harmonic potential. In the bubble state the wave 
function is localized around the surface of a sphere so 
that the matter density forms a spherical shell or bub- 
ble. We give a detailed account of the preparation of 
bubble states and their decay rates induced by nonadia- 
batic leakage. In particular, we address the question of 
how to compensate for gravitational effects in the labora- 
tory which otherwise would impede the creation process. 
We expect matter- wave bubbles to have interesting appli- 
cations some of which are investigated in Sec. IV. There 
it is shown that they can be used as stepping-stones to 
prepare atoms in highly excited eigenstates of the bare 
harmonic trapping potential and to create new types of 
non-linear eigenstates of BECs. Furthermore, we exam- 
ine collapse and revival effects which are found when the 
external coupling is switched off. By also switching off 
the trapping potential we can observe the free bubble 
expansion. 

In Sec. V we turn to the investigation of gravitational 
influences on the bubble preparation process. At first 
sight gravity is seen to be detrimental, but it is shown 
that it can be exploited to obtain a novel scheme for the 
creation of two-dimensional BECs. Gravity will cause 
the atoms in the shell potential to pool at the bottom 
of the trap so that the condensate forms a flat disk. In 
Sec. V the initial transfer of a BEC into such a disk state 
in the course of switching on the external field is exam- 
ined numerically. In particular, we show that if the field 
parameters are changed appropriately the disk radius is 
steadily increased whereas its width keeps shrinking due 
to enhanced confinement. This suggests that the con- 
tinuation of this process could ultimately yield a two- 
dimensional trap for the BEC, although the numerical 
study of the approach to this limit is not currently feasi- 
ble. In Sec. VI we further pursue this idea by giving some 
general qualitative, and semi-quantitative, estimates and 
arguments about the conditions necessary to reach the 
2D regime. They indicate the feasibility of our approach 
by showing that the requirements for the applied fields, 
preservation time of the condensate etc. are demanding 
but still within the reach of currently available technol- 
ogy. The paper ends with brief conclusions and outlook 
given in Sec. VII. 

The most important difference to the previous propos- 
als of Refs. [7, 8] for creating 2D traps is the fact that our 
method relies on adiabatically deforming a conventional 
magnetic trap and does not require incoherent processes, 
e.g., optical pumping, for loading. This would allow for 
working with extremely cold, dense, and, possibly, coher- 



ent atomic ensembles throughout the whole process. In 
comparison to the experiments [4, 5, 6], our scheme does 
not make use of optical potentials, but a combination of 
magnetic and RF fields. It thus avoids potential diffi- 
culties with spontaneous emission in very steep optical 
traps that require high laser intensities. In contrast to 
Refs. [5, 6], our proposal produces a single condensate 
with a large number of atoms, similar to the experiment 
of Ref . [4]. However, the trapping frequencies obtainable 
in the RF-scheme might be significantly higher trapping 
frequencies than the ones reported in the latter work. 



II. ADIABATIC POTENTIALS 

The basic principles of adiabatic potentials can be 
understood from examining a quantum-mechanical two- 
state particle that propagates in the vicinity of a linear 
potential crossing. In the interaction picture, with re- 
spect to the coupling field that gives rise to the crossing, 
the Schrodinger equation for this system is written as 
follows: 
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We can transform this equation to a basis that diago- 
nalizes the potentials ±Cr and the coupling Q at each 
point, i.e., to the basis of the spatially dependent 'dressed 
eigenstates.' In this basis the Schrodinger equation has 
the form [11] 
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The potentials V± (r) are the adiabatic potentials and are 
given by 



V±(r,t) =±^{Cr) 2 +fl 2 . 



(3) 



They arise from the pointwise diagonalization of the 2x2- 
matrix formed by ±CV and J7. The terms 
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stem from the nonlocal character of the kinetic energy 
term. 

We see that V+(r) + Vfej n (r) is a binding potential. 
If the coupling term V c (r, d r ) were not present, a wave 
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packet prepared in the state '+' would remain trapped 
forever in the potential V + (r) + Vki n (r). However, if £1 
is small the coupling V c is dominant. In this case any 
wave packet quickly leaves the crossing region around 
r ~ 0, and the description given by Eq. (2) is not useful. 
Nevertheless, if fi is increased, the kinetic coupling V c (as 
well as Vkin) rapidly becomes small [24] and the motion 
of the wave packet is more and more determined by the 
adiabatic potentials V±. In fact, it is shown in Sec. III.B 
that the lifetime of a wave packet prepared in the internal 
state '+', i.e., the time it takes the wave packet to transfer 
to the state '— ' and leave the crossing region, increases 
exponentially with Q. 

At first sight, it seems counterintuitive that it is possi- 
ble to trap a particle, even with strong coupling, between 
two potentials in a region where the particle is not sta- 
ble. We can draw an analogy with a particle moving in a 
magnetic quadrupole field. Each of its bare Zeeman sub- 
states (with bare meaning having fixed spatial orienta- 
tion) individually would experience the field as unstable. 
However, if the particle moves slowly enough, couplings 
between the states are induced so that the particle's ori- 
entation with respect to the local direction of the total 
magnetic field is preserved. In this way, trapping can en- 
sue for weak-field seeking states. In the case considered 
here, a slow particle tends to remain in the same dressed 
eigenstate and its motion is governed by the adiabatic 
potentials. 



III. CREATION OF MATTER- WAVE BUBBLES 
A. Basic approach 

To work out the basic ideas of our approach we first 
discuss its realization in the absence of gravity. As a 
result, we obtain a scheme to produce thin, highly stable 
matter- wave bubbles or shells in which the trapped atoms 
are localized around the surface of a sphere. 

The starting point for our method is a coherent sample 
of atoms produced, e.g., by Bose- Einstein condensation 
and trapped in the ground state of a harmonic magnetic 
potential. The preparation scheme then proceeds by ap- 
plying a sequence of radio frequency fields that couple the 
initial internal atomic state, a weak-field seeking Zeeman 
sublevel, to a second hyperfine ground state (multi-level 
excitation schemes can be considered as well, see Sec. 
V). One interesting aspect of the preparation process is 
that it may appear to be counterintuitive as the second 
state is an untrapped, high-field seeking state [see in- 
set, Fig. 1(a)]. As this technique also forms the basis of 
evaporative cooling, one may be led to expect that this 
procedure will inevitably cause a rapid depletion of the 
trapped atomic population. Nevertheless, as was antici- 
pated in Sec. II and is also shown below, if the fields are 
controlled in an appropriate way one may also obtain a 
very different effect. 

To model the creation of our bubble states we study 
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FIG. 1: (a) Schematic of field- induced adiabatic potentials V± 
for A > 0. Dashed curves show the bare potentials crossing 
at r c . Inset: bare potentials showing resonance at r c . (b) 
Bubble, or shell state as obtained by the preparation scheme 
discussed in the text with (in scaled units) Af ina i = 60.0 
and £1 fi na i = 9.0 (see inset). Full curve: atomic density 
|0+| 2 /r 2 in the adiabatic state +; dotted and dashed curves 
show |0i ] 2 /r 2 and \<p2\ 2 /r 2 , respectively. Nonlinear effects are 
not included here. 

the coherent time evolution of a condensate initially pre- 
pared in a hyperfine state |1) in the ground state of a 
spherically symmetric magnetic trap. When t > ex- 
ternal fields are applied that induce a coupling of state 
|1) to a Zeeman sublevel |2) whose magnetic moment is 
supposed to be equal in magnitude but opposite in sign. 
Such a coupling may be realized, e.g., as a transition be- 
tween the hyperfine ground states \F = l,mp = —I) 
and \F = 2,m F = -1) in 87 Rb [12]. The field-induced 
coupling strength is denoted as O and is spatially inde- 
pendent to a good degree of approximation. In the fol- 
lowing model we assume a radial I = form of the wave 
function components ipi{v) — <pi{r) / t \[knr : Then work- 
ing in an interaction picture with respect to the applied 
fields, the time development is determined by the radial 
Gross-Pitaevskii equation 

. / 1 d 2 r 2 A(i)\ 1 

N (Un\fa\ 2 + U 12 \<f> 2 \ 2 ) ^, 
. / 1 d 2 r 2 A(t)\ 1 

N{U l2 \^\ 2 + U 22 \<j> 2 \ 2 )^. (6) 

In this equation time is scaled to the trap inverse os- 
cillator frequency length is scaled to the har- 
monic oscillator length scale a ho , and fl(t) and A(t) 
are scaled to Tilo. The effective detuning is defined by 
A(t) = [hcjf - AE(0)]/Jkj, where AE(0) is the en- 
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ergy difference between the two hyperfine states at the 
origin (trap centre) and ujf the frequency of the ap- 
plied field. The nonlincarity parameters are given by 
Uij = dij/aho with ciij the scattering lengths for intra- 
and inter-species collisions. The states are normalized 
according to J °° dr(\4>\\ 2 + \4>2\ 2 ) = 1- The total number 
of atoms is denoted N. 

The strategy that we pursue in our engineering scheme 
is to control the condensate by slowly changing field- 
induced adiabatic (or dressed) potentials. These poten- 
tials, which are defined as the spatially dependent eigen- 
values of the potentials and couplings in Eqs. (6), are 
given by [cf. Eq. (3)] 

V±(r,t) =± V /P r= A(i)] 2 /4 + fl 2 (t) (7) 

and are depicted for A > in Fig. 1(a). The potential 
V- actually gives rise to the evaporative cooling effect in 
the usual arrangement. In that case one applies a field 
with an effective detuning which is large compared to the 
mean particle energy. The atoms then move in the po- 
tential V_, and the ones that reach its maximum, from 
the left, with sufficiently slow velocity, go over the top 
and get expelled from the trap. However, we will show 
that the atoms can also be prepared in the lowest-energy 
eigenstate (or, more exactly, resonance) |0) = |0; O, A) of 
the potential V + . This quasi-bound or 'trapping' state 
|0) will realize the spherical shell state or matter- wave 
bubble (see Fig. 1(b) and Rcf. [13]). The state will be lo- 
calized around the crossing of the two bare potentials at 
r c = \/A and have a width of Ar = (fi/A) 1 / 4 , provided 
a harmonic expansion around the potential minimum is 
justified. The state is a genuine superposition of the in- 
ternal states |1) and |2). 

B. Lifetime of bubbles 

Before considering the bubble preparation process in 
detail one question immediately arises, i.e., the stability 
of the system once it is prepared in the state |0). At first 
it is not obvious that atoms may remain trapped for a 
substantial time at the point of maximum effective cou- 
pling between states |1) and |2). As we have discussed in 
Sec. II, it becomes more plausible if one transforms Eqs. 
(6) to the dressed state basis, i.e., the basis that diagonal- 
izes the bare potentials and the coupling at each point r. 
In this picture the two wave function components appear 
coupled by kinetic terms whose significance is rapidly 
diminished when is increased. If non-linear interac- 
tions can be neglected, the decay rate 7 of the trapping 
state |0) can be determined with the help of semi-classical 
methods developed in connection with molecular prcdis- 
sociation [14, 15]. As we show in the Appendix, from 
these techniques it follows that 7 = — 2Im£'o, where the 
complex ground state energy E is determined as a solu- 
tion of 

^S(E) _ l] cos $(£) e -i[/3(£)-*(i?)] + cos /?(£?) = (8) 
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FIG. 2: Evaluation of decay rates 7(fi) for bubble states 
according to Eqs. (8) (full curve) and (9) (dashed) at A = 
61.0, using the approximations of Ref. [14]. Circles: numerical 
values determined from Eq. (6). 



with /3(E) = tt(2E + A - l)/4 and the parameters 6(E) 
and Q(E) characterizing the scattering matrix of the lin- 
earized potential crossing problem. For these quantities 
there are several analytical approximations in the liter- 
ature [16]; following, e.g., Ref. [14] one can put 6(E) = 
l/8aband<f>(E) = 2b 3 /3a + avgT[i6(E)]+6(E) \n[6(E)}- 
26(E)\u(b/a) + tt/4 with a 2 = A/(8ft 3 ) and b 2 = E/Sl. 
For large enough one obtains 

2cos2 ^) (9) 
' {cx P [2tt6(E)] - l}(d$/dE) 1 ' 

where all qu antities have to be evaluated at ReE — 
fi + \J A/451. The comparison in Fig. 2 between the 
predictions of Eqs. (8) and (9) and the direct numeri- 
cal determination of decay rates from Eqs. (6) indicates 
the validity of these approximations. To obtain the nu- 
merical decay rates we first generated the bubble state at 
a required (A, SV) by performing the preparation process 
described in Sec. III.C. We then monitored the decay of 
the norm of the wavefunction in the internal state '+' as 
a function of time. After an initial transient, this decay 
was exponential to a very good degree of approximation. 

Equation (9) yields two important insights: firstly, the 
decay is exponentially suppressed with growing f2. In 
the limit of 7 <C 1, this exponential suppression may be 
approximated as 

lex ~ 2 exp(-^r! 3 / 2 /V2A 1 / 2 ). (10) 

Secondly, for ReE = 2k + 3/2- A/2, with integer k, the 
decay rates become very small. In these cases the state 
|0) is in resonance with an eigenstate of the bare har- 
monic trapping potential. This stabilization effect may 
be used to obtain extremely long-lived states already for 
moderate coupling strengths. 
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C. Preparation 

Having established the existence of long-lived spheri- 
cal bubble states in adiabatic potentials, we now turn to 
the details of their preparation. The first phase is to let 
the bare harmonic potential of state 1 evolve smoothly 
into the dressed potential V + and the second phase is to 
expand the bubble outwards. Both phases are done so 
that the wave function is kept at all times in the instanta- 
neous ground state |0; O(t), A(t)). In the example of Fig. 
1(b) (in which non- linear effects have been neglected) 
we achieve the first phase by increasing to the final 
(scaled) value of 9 at fixed negative detuning A = — 3 
(see Fig. 1(b), inset). Then the second phase of bubble 
expansion is achieved by maintaining the intensity, but 
gradually ramping up the detuning. In this way less than 
4% of the initial population was lost before reaching the 
final state shown in Fig. 1(b). The particular merit of 
this approach is that there is no need to follow a pre- 
cisely prescribed path as the method relies on adiabatic 
following (in a similar way to the wave packet guidance 
in APLIP [17]). Any field sequence that guides the sys- 
tem sufficiently slowly through states of long lifetime is 
suitable. 

If we make our atomic bubbles from an initial BEC, we 
should also study the effects of the non-linear interactions 
introduced in Eq. (6). The dressed-state description is 
readily generalized if Uu ~ U22 — Uu — U (such a 
situation is realized, e.g., in the Na F = 1 hypcrfinc 
multiplet [1]). Under these circumstances we can still 
define an adiabatic wave function <f> + that evolves, to a 
good degree of approximation, according to the Gross- 
Pitaevskii equation 



-l^2+V + (r,t) + NU 



\4>+ 



|21 



(11) 



In particular, for a sufficiently strong nonlinearity U, the 
atomic density can be described by the Thomas-Fermi 
ansatz 



\H-V+{r,t)]/NU 



(12) 



with the chemical potential fi determined by the nor- 
malization condition. Numerical studies show that shell, 
or bubble, states can then be manufactured in the same 
way as before. Once prepared, their decay is now non- 
exponential, but the rate is of the same order of magni- 
tude as the corresponding system in the linear case. If 
the atomic cross- and self-interactions differ significantly, 
numerical simulations indicate that the creation of bub- 
ble states becomes more difficult and their lifetime is re- 
duced. 



will cause the atoms to pool at the bottom of the shell 
potential. In this subsection, we first give a simple esti- 
mate of how strong the gravitational influence must be 
to cause significant deviations from a homogeneous den- 
sity distribution in the bubble state. We then propose a 
method to compensate for gravitation in the laboratory 
with the help of optical potentials. 

When the bubble width is small compared to the ra- 
dius, the radial dependence of the adiabatic potential is 
approximately harmonic with a minimum at the bubble 
radius r = VA and an angular frequency luo = \J A/Q, 
[13]. Thus, we can approximate the Hamiltonian for 
atoms in the trapping internal state '+' as 

V 2 1 

H = -— + -Lu (r-ra) 2 + Gr cos 6, (13) 

with the radial and polar coordinates r, 0, an d the scaled 
gravitational acceleration G = g\Jmj hiof (with g ~ 
9.81ms -2 ). We now make the ansatz 

^(^)^^ CXPh(r ' r ° )2/2fT2] /(cos^) (14) 



for the ground state of this Hamiltonian, where a = 
1/y/uJo. This ansatz implies that, radially, the wave 
function is always in the ground state, i.e. the influ- 
ence of gravity is manifest only in the polar envelope 
/(cos#). For present purposes it is sufficient to choose a 
very simple form for /, e.g., /(cos#) = acos8 + b. The 
normalization factor in Eq. (14) is then given by Af — 
y/n(2a 2 /3+b 2 )/a. Furthermore, a variational calculation 
shows that a and b are related by a/6 = Q — w 'Q 2 + 3/2 
with Q = 3/4GrQ. A significant influence of gravity is 
certainly present when /(cosO) = 0, i.e., a/b = —1, as in 
this case the bubble has 'opened up.' This is realized for 
Q ~ 0.17 or, in unsealed coordinates, 



r ~ (3^/gm 2 ) 1 / 3 



(15) 



which is about 5.5 x 10~ 7 m for Rb87. In the presence of 
gravity, bubbles can therefore only be observed in spher- 
ical traps with trapping frequencies large compared to 
400Hz. The estimate (15) can be interpreted as a bal- 
ance condition between kinetic and gravitational energy 
which are of the order h 2 /2mr 2 and mgro, respectively. 

The detrimental influence of gravity can, at least in 
principle, be compensated for by exposing the trapped 
atoms to an additional optical dipole potential. In a typ- 
ical Gaussian beam configuration this potential, which 
acts on all hyperfine sublevels in the same way, is given 

by [1] 



V d (x, g) = 



1 + (x/x R y 



■ cxp 



2g 2 



w 2 [i + ( x /x R y 



(16) 



D. Bubbles with gravitational field compensation with Q = y / y 2 + z 2_ The R ay l e igh length x R and the 

beam waist radius Wo are related by x R — ttwq / X op with 
So far, we have neglected the influence of gravity. Its \ op the wavelength of the optical field. Finally, in un- 
effect on the bubbles is expected to be detrimental as it scaled units, Vo = 3Tc 2 P/ AwqOJq where ujq denotes the 
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that in this direction the variation of the combined op- 
tical and gravitational potential is almost negligible. To 
find an optimized shape for given f2 and A, it is advis- 
able to experiment by slightly varying the laser power 
and the position of the magnetic trap center. Neverthe- 
less, the preparation scheme is robust in the sense that 
(for the configuration described above) one still obtains 
a bubble if the location of the magnetic trap (in z) rel- 
ative to the dipole potential is changed by about ±10%, 
and the laser intensity by about ±2.5%. However, during 
preparation these parameters have to be stabilized very 
well (to within a fraction of the indicated range) to avoid 
excitation of the bubble. 



FIG. 3: Schematic showing the local compensation of gravity 
with the help of an optical potential. Shown are dipole po- 
tential Vd, the gravitational potential V g , and their sum. The 
arrow indicates the flat region of the combined potential in 
which a bubble might be created. 



resonance frequency, T the spontaneous decay rate of 
the excited state, A the detuning between laser and the 
atomic transition and P the applied laser power. The 
idea is now to adjust the trap depth Vq such that at 
the turning point z t of the optical potential along the 
line x = y = 0, which is defined by d 2 Vd(z t )/dz 2 = 0, 
the slope dVd(z t )/dz = —mg. In this way, the com- 
bined optical and gravitational potential is almost con- 
stant around z t , the lowest-order corrections being cubic 
in z — z t (see Fig. 3). The magnetic trap has then to 
be placed inside this area. From Eq. (16) it follows that, 
as z t — wo/2, we have to choose Vb = exp(l/2)w TOg/2 
to fulfill the above condition. If the condensate is in 
the Thomas-Fermi regime, the z extension of the volume 
within which a bubble could be created is somewhat less 
than (wQn/mg) 1 / 3 with fi the BEC chemical potential 
[for which an estimate is given in Eq. (20)]. This follows 
from stipulating that the variation of the combined op- 
tical and gravitational potential over this volume should 
be less than fi. 

We have performed numerical simulations to demon- 
strate that this approach actually offers the possibility to 
overcome the influence of gravity. To this end, we have 
studied the preparation process using a two-dimensional 
version of Eq. (24) including the optical potential. Equa- 
tion (24) describes time evolution in the effective adi- 
abatic potential of the Rb87 F — 2 multiplet. As an 
example we have shown that in an anisotropic trap with 
v z — w z /2-7r = 220Hz, v x = 30Hz, bubble states for a 
BEC corresponding to an atom number of 10 5 in three 
dimensions can be produced at A = 6.6kHz, ft = 2.6kHz, 
Vb = —129kHz, and w = TSfim. For the preparation, ft 
is switched on within At = 29ms at A = — 1.1kHz, and 
then A is increased to its final value within At = 200ms. 
The bubble has radial and axial diameters of 15^m and 
110/im, respectively. The necessary laser power is esti- 
mated at around 1W. The large axial extension shows 



E. Bubbles and trap anisotropies 

At this point it should be mentioned that — at least in 
the absence of non-linear interactions — the preparation 
of bubbles is impeded by another problem besides gravity, 
namely, the unavoidable trap anisotropies. Imagine that 
in the bare magnetic trap u> x — (3u> z , but we keep /3«1. 
Then the effective trapping frequency along the x axis is 
given by (i^J A/9l (scaling now with respect to u) z ). We 
can model our system Hamiltonian as 



(17) 



where u>o(0) appropriately interpolates between \J A/Q, 
and 0y/ A/£l. The ground state is approximated as 



'tr 



1 expKr-r ) 2 /2a(g) 2 ] 



/(cos0) (18) 



with a(6) = 1/V^oW We again could determine 
/(cos 9) by minimizing the energy functional. However, 
following the discussion of the previous paragraph we can 
argue that variations in the angular kinetic energy carry 
a cost of the order of 1/ A whereas the difference in 'polar 
potential energy' is given by \/3 2 — l|A/f2. Therefore, we 
expect trap anisotropies to have a significant influence as 
soon as 



A 
9. 



> 



1 



(19) 



As typically A ^> 1 and A/Q 3> 1, the trap frequencies 
have to be adjusted very carefully in order to enable the 
generation of bubbles. In the case of Fig. 2, for example, 
the estimate (19) indicates that \uj x /lu z — 1| should be 
of the order of 0.1%. Fortunately, the problem can be 
alleviated to some extent by making use of the atomic 
interactions. In the Thomas- Fermi limit, the chemical 
potential of a thin bubble is estimated to be 



/ 3 UN \ 



2/3 



(20) 
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with U 



it/o-ho as the scaled nonlinearity coefficient 



and where TV is the number of atoms in the bubble. 
Thereby, /x is measured with respect to the bottom of 
V+ (r) . We expect that only if the non-linear interaction 
energy — which is of the order of /x — is in sufficient excess 
of the effective potential energy, i.e., if 



li 3 / 2 - 



(21) 



then a bubble of approximately constant density can be 
formed. Thus, even if gravity is compensated for, bubbles 
can probably only be produced in the presence of strong 
non-linear interactions. On the other hand, the numer- 
ical simulations of the previous subsection have shown 
that under such conditions a preparation may indeed be 
possible. 
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IV. APPLICATIONS OF BUBBLES 
A. Production of excited harmonic oscillator states 

As a first application of matter-wave bubbles we now 
discuss the preparation of excited harmonic oscillator 
states in the absence of non-linear interactions (i.e., for 
sufficiently dilute samples). If, after creating a bubble 
state, the coupling strength is slowly reduced at fixed A 
the system will again evolve through a sequence of instan- 
taneous eigenstates. This time, however, the eigenstate 
will not significantly change its energy relative to the 
minimum of the bare trapping potential. Qualitatively 
speaking, the wave function gradually moves out from 
the crossing region where it is being trapped and begins 
to experience more strongly the presence of the bare har- 
monic potential. This process is accompanied by a loss of 
atoms which end up on the repulsive potential after leav- 
ing the crossing region. At £1 = the remaining bound 
wave function will have reached an excited harmonic os- 
cillator eigenstate which, in the spherical harmonic case, 
is characterized by having I — 0. As mentioned before, 
the energy of this eigenstate is approximately equal to 
the energy of the initial bubble state. To give an exam- 
ple, if we start from the state shown in Fig. 1(b), and the 
coupling strength is ramped down within a time interval 
of At = 16, we arrive at an eigenstate of energy 35.5 
with respect to the minimum of the harmonic potential 
[13]. A sequence of intermediate states appearing in the 
course of this process is shown in Fig. 4. After comple- 
tion, the admixture of other eigenstates is less than 2%, 
and the population is 36% of the initial population in the 
harmonic trap ground state. 

A few remarks should be made about this preparation 
scheme: (i) It has to be performed 'quasi-adiabatically', 
i.e., slow enough to avoid substantial excitation of other 
eigenstates, but sufficiently fast to reduce losses as much 
as possible, (ii) Although the system effectively ends up 
with the same Hamiltonian as in the very beginning (cou- 
pling switched off) and the time evolution is performed 



FIG. 4: Resonance states of Eqs. (6) at A = 60 and various 
values of £1 (nonlinearities not taken into account). The wave 
functions are determined numerically by slowly decreasing the 
value of £1 in Eqs. (6) after initially preparing a bubble state. 
Full curves: |0i| 2 , dashed: <j>2 1 2 , dotted: sum of both. The 
displayed wave functions are normalized to one. 

quasi-adiabatically, the final state is very different from 
the initial one. This is because the Hamiltonian depends 
on two external parameters Jl and A, and O = rep- 
resents a singularity in the sense that the value of A 
becomes irrelevant, (iii) The oscillator eigenstates can in 
turn be used as intermediate states to produce radially 
excited bubble states. To this end, one simply switches 
the field coupling on again with a reduced detuning. 

B. Nonlinear eigenstates 

The bubble states can be regarded as the ground states 
of a specifically tailored potential. Recently, however, the 
study of macroscopically excited states of BECs has re- 
ceived much attention [18, 19] and the manufacture of 
harmonic oscillator eigenstates, as presented above, indi- 
cates a way to prepare a new class of such excited states. 
This class can be thought of as the non-linear general- 
ization of the eigenstates of linear systems; for a one- 
dimensional model some properties of such 'non-linear 
modes' were examined in Ref. [20] without discussing 
ways for their actual preparation. For highly excited 
states, non-linear effects are expected to play a minor 
role, in general, because of the reduced density. We thus 
focus here on the first excited non-linear I = mode 
which is characterized by one radial node in the wave 
function. In the absence of non-linear interactions the 
preparation scheme of Sec. IV. A works equally well for 
low- and high-lying eigenstates, so we have numerically 
applied the same approach, with suitably chosen values 
for £1 and A, to the full Gross-Pitaevskii equation. Our 
studies show that the scheme is still applicable, though 
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FIG. 5: Condensate density \(f>i\ 2 /r 2 at times At = 2, 12, and 
32 after completion of the preparation scheme. Preparation 
proceeds from the bubble state at Q, = 5.6 and A = 5.4 
by reducing Q to within At = 70. The bold curve shows 
the exact 'non-linear mode' for the final value NUu = 17.2. 
The inset compares this state (dashed) to the corresponding 
eigenstate in the absence of nonlinearities (full curve). 



the atomic interactions cause detrimental effects: shortly 
before reaching the final value of f2 = one encounters a 
very strong loss of atoms, and, because of this, the pro- 
cess has to be performed rather quickly resulting in an 
appreciable excitation of the final wave function. Fur- 
thermore, the final detuning A has to be selected more 
carefully (to within ±0.1 u) to optimize the number of 
atoms that remain trapped. In spite of these difficulties, 
satisfactory results can be obtained for final nonlinearity 
parameters NUu up to the order of 10. As an exam- 
ple, Fig. 5 shows the density distribution \<f>i(r)\ 2 /r 2 at 
various instances after completion of the preparation pro- 
cess (for which J7 n = U 2 2 = I.O8U12 was assumed, as in 
23 Na). The breathing seen in the wave function indicates 
that there is additional radial excitation. The parameter 
NUu = 17.2, and the projection onto the exact station- 
ary state at this value varies between 60% and 95%. The 
efficiency, i.e., the ratio of final and initial atom number, 
is 9.5%. The comparison with the eigenstate of the linear 
case (see inset) shows a broadening of the wave function 
and a significant reduction of the central density due to 
the interatomic repulsion. 



C. Collapse, revival, and free expansion 



FIG. 6: Condensate density \<f>i\ at times At = 1.6 and 3.2 
(inset) after suddenly switching off the RF coupling Q. The 
different curves have nonlinearity parameters UN — (full 
curve), 5 (dashed), and 20 (dotted). In each case, the initial 
state was the respective matter-wave bubble at A = 60 and 
Q = 9. 



this scenario repeats itself periodically, the wave function 
regaining its initial shape at times mr, n — 1,2,... In 
the presence of atomic interactions, however, the shape 
is gradually distorted. Figure 6 shows |</>i(r)| 2 for three 
different values of gN at times t = 1.6 and 3.2, i.e., for 
'complete' collapse and re-expansion, respectively. In the 
large figure the wave functions show an interference pat- 
tern which is due to particles from opposite sides of the 
bubble passing through each other. Note that for grow- 
ing non-linear interaction the central interference fringes 
are pushed outwards in agreement with earlier studies on 
similar systems (see, e.g., [21]). It might be possible to 
infer the nonlinearity parameter from the fringe pattern. 
The inset in Fig. 6 illustrates the subsequent broadening 
of the wave function in the course of the re-expansion 
due to non- linear interactions. 

A related behavior can be observed if the magnetic 
fields are switched off along with the RF coupling. As it 
is no longer subject to any potential, the localized radial 
bubble wave function displays a time evolution similar 
to a free particle and gradually broadens due to disper- 
sion (see Fig. 7). After some time wave function pieces 
from opposite sides of the bubble start to overlap each 
other and an interference structure ensues. Again, the 
interference fringes are shifted outwards if non-linear in- 
teractions become significant (see the inset). 



Another interesting effect occurs if, after creating the 
bubble, the coupling strength is instantaneously re- 
duced to zero. In this case, the two components of the 
bubble evolve almost independently of each other in their 
respective bare potentials. Component 2 is therefore 
rapidly expelled from the trapping region, whereas com- 
ponent 1 undergoes a collapse or contraction into the 
center of its binding harmonic potential followed by a 
re-expansion. If non-linear interactions can be neglected 



V. DISC-SHAPED CONDENSATES 

In this section we examine how our trapping scheme 
is modified in the presence of gravity and in the absence 
of any gravitational compensation (such as described in 
section HID). The results will also lay the foundation 
for our proposal for the creation of 2D atom traps which 
is outlined in the following section VI. 
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FIG. 7: Condensate density |0i| 2 + \4> 2 \ 2 at times At = 0.0 
(long-dashed), 2.0 (full), 4.0(dotted), and 6.0 (short-dashed) 
after suddenly switching off both the RF coupling and the 
magnetic potential. The initial state is the matter-wave bub- 
ble at A = 40, O, = 9, and UN — 0. In the inset the conden- 
sate density is shown at At = 6.0 for nonlinearity parameters 
UN = (full curve), 20 (dashed), and 40 (dotted). 



A. Adiabatic trapping in the presence of gravity 

To make the presentation more concrete, we look at a 
specific example that corresponds to typical experimental 
conditions. We thus consider a condensate of about 10 5 
Rb87 atoms initially prepared in the F = 2, Mp = 2 hy- 
pcrfine sublevel and in the ground state of an anisotropic 
magnetic trap with frequencies v y = v z = 220 Hz and 
v x = 11 Hz. The z axis again points along the vertical 
direction. An RF field is then applied that couples the 
sublevels within the F = 2 multiplct. By appropriately 
tailoring the time dependence of the field, the BEC al- 
ways remains in the ground state of the RF-induced adi- 
abatic potential and is thus manipulated in a controlled 
way. 

The condensate dynamics is determined by the Gross- 
Pitaevskii equation for the components ipM, M = 
2, . . . , 2, 



iipM = 



V 2 M 

-— + —{k 2 x 2 + y 2 + z 2 ) - MA(t) + Gz Vm 

+^M,M-l(t)V'M-l + ^M+l,M(t)lpM+l 

+UNp(r,t)^ M - (22) 



Again, all quantities arc dimensionless. They are now 
scaled to units derived from the radial angular frequency 
uo z = 2-kv z . The effective detuning is given by A(t) = 
[hujf — AE(0)]/hio z with AE(0) the Zeeman energy split 
between subsequent hyperfine sublevels at the magnetic 
field minimum. The parameter k — v x jv z gives the ra- 
tio between axial and radial trapping frequencies. The 
scaled gravitational acceleration G becomes 7.06 in the 
present setup. The R F co upling constants are given 
by n 2 ,i = -2 = y/2/3Q lfi = y/5/3tlo,-i = Sl(t). 
For simplicity, all non-linear interactions coefficients are 
taken to be equal and denoted by U; N is the total num- 



As discussed in Sec. III.C, the simplified form of 
the non-linear interaction allows us to transform to the 
dressed-eigenstate basis in the same way as in the lin- 
ear case. The field-induced adiabatic potentials are now 
given by 



V M (r, t) = M^/[V tr (r)-A(t)] 2 + W(t), 



(23) 



where V tr {r) = (k 2 x 2 + y 2 + z 2 )/4, so that a wave packet 
initially prepared in the state M = 2 can evolve in the 
potential V 2 (r, t) + Gz if the potential is deformed slowly 
enough. The approximate equation of motion for ip 2 thus 
reads 



ii> 2 = 



-— + V 2 (r,t) + Gz + UN\i> 2 (r,t)\ 2 



$24) 



ber of atoms and p(r, t) = 



M = 



\^M{r,t)Y 



The influence of gravity becomes apparent by writing 
the z-dependent part of the potential as Mz 2 /4 + Gz — 
M(z + 2G/M) 2 /4 - G 2 /M. We see that, effectively, 
the position of the bare trap minimum is shifted to 
z = —2G/M (which typically is large compared to the 
ground state extension), the minimum itself is shifted by 
-G 2 /M. At A = all bare potentials touch at z = 0. 
The bare potential for M — —2 intersects the M — 2 po- 
tential at its minimum for A = G 2 /4. As the condensate 
is initially localized at this minimum, we conclude that 
only for A > G 2 /4 does a significant shift and deforma- 
tion of the ground state set in. For large A the ground 
state is located at about 2y/A. In Fig. 8 we show cuts 
along the potentials in the z direction that illustrate the 
characteristic behavior explained above. 

What does the full spatial dependence of the poten- 
tial 1^2 (r, t) + Gz look like? The potential V 2 (r,t) alone 
has its minimum on an ellipsoidal surface defined by 
Vt r (r) = A(t). However, due to the influence of grav- 
ity the combined potential V 2 (r, t) + Gz is strongly tilted 
in the z direction, so that the wave function will assemble 
around its bottom (we will see this in Fig. 10). The main 
control parameter to vary the shape of the adiabatic po- 
tential is the detuning A. At a given A the coupling 
strength SI has to be chosen large enough that the life- 
time of the condensate becomes sufficiently long. 



B. Preparation 

To demonstrate that atoms can be actually be trapped 
in the RF-induced adiabatic potential we have performed 
numerical simulations of Eq. (22) in two dimensions, i.e., 
the direction of gravity z and the weak trapping direction 
x. These calculations should be able to capture the main 
aspects of the wave packet behavior. Our results indicate 
that the same two-step approach as outlined in Sec. Ill 
can be used to transfer the BEC into the adiabatic trap. 
As an illustration, we show in Fig. 9 the result of one of 
our simulations. For the 2D calculation, the nonlinear- 
ity parameter NU was chosen so that the extension of 
the ground state in the bare magnetic trap [which is the 
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Bare trap potentials 
without gravity 



Bare trap potentials 
and gravity mgz 
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FIG. 8: Schematic of the potentials along the z-direction 
showing the characteristic influences of gravity and the ap- 
plied RF-field. 



initial state ip2{t = 0)] coincides with the extension in x 
and z of the 3D ground state for the given atom number 
of 10 5 . In the first step of the preparation scheme, the 
RF intensity is linearly ramped up to the desired final 
value of ft at negative A. In Fig. 9, in order to reach 
a final CI = 12 (rs 2.64kHz) the field is switched on at 
A = -5(w -1.1kHz) within At = 20 ~ 14.4ms. In the 
second step the RF detuning is simply increased to the 
final value thereby keeping the intensity fixed. In Fig. 9, 
the detuning A was increased to 60(~ 13.2kHz) within 
At = 160 ~ 116ms. The simulations already show the 
deformation effect that we build upon in the next sec- 
tion to obtain 2D trapping, namely, a squeezing in the 
z direction indicating tighter trapping, and an expan- 
sion in x showing the decrease in the corresponding trap 
frequency. Furthermore, we see that the Thomas-Fermi 
approximation gives a very good description of the wave 
function. Nevertheless, it should be mentioned that dur- 
ing preparation a slight excitation of the wave function 
along the x axis appears (not visible in Fig. 9) indicating 
some nonadiabatic effects. 

The effort required for the numerical calculations 
shows that it is very difficult to go significantly beyond 
the trap deformations shown in the above example — and 
in particular to simulate reaching the 2D limit. In the 
next section we give some semi-quantitative arguments 
to determine the values of CI and At that should be used 
to reach a given A. Here we restrict ourselves to a qual- 
itative discussion. The chosen A places a lower limit 
on the value of CI — only for a high enough CI is the life- 
time of the BEC in the adiabatic potential sufficiently 
long. For this complicated five-state system it is difficult 
to estimate the lifetimes analytically in an accurate way. 
However, the numerical calculations indicate that with, 
for example, A = 60, a value of CI = 12 should be ad- 
equate. Furthermore, the study of the linear two-state 
system shows that the lifetime increases exponentially 
with CI. We thus expect that the lifetime can always be 
adjusted by a moderate increase in Cl. 

Concerning the rise times At, they are bounded from 
below by the requirement that we want to avoid nonadi- 




FIG. 9: Loading an adiabatic trap in the presence of 
gravity. Shown is the ^-integrated density P(z) = 
J dx X^m-_2 ( x > z )| 2 a t the indicated times for the prepa- 
ration process described in the text. The inset shows the z- 
integrated density P(x) at t = and 180. The dashed curves 
show P as obtained from the Thomas-Fermi approximation. 



abatic excitations and keep the BEC in the ground state. 
We expect that it is more likely we would generate such 
excitations along the x direction than along the radial di- 
rections as in the former the eigenstates are spaced much 
more closely. In fact, this is observed in the numerical 
simulations as mentioned above. The amount of excita- 
tion is reduced by performing the process more slowly. 
An upper limit on the rise time is placed by the lifetime 
of the adiabatic ground state and the overall decay time 
of the BEC. 



C. Signatures of adiabatically trapped BECs 

There are various signatures which allow one to verify 
experimentally that the condensate behavior is indeed 
determined by the adiabatic potential (apart from the 
fact that the BEC is still there at all): 

1. The vertical position of the condensate is changed 
(lowered). With z = indicating the minimum of the 
magnetic trap field the condensate is located around 
z = — G (~ — 5.1/xm for our example) in the absence 
of the applied RF fields. With these fields turned on at 
sufficiently large A the potential minimum is shifted to 
the vicinity of the intersection between the different bare 
potentials at x = y = 0, z = — 2\/A. For A = 60, e.g., 
the condensate center is transferred to — 11.3/um, i.e., it 
is shifted over a distance larger than the BEC extension 
in the z direction. 

2. The shape of the condensate is changed. At A = 40, 
Cl = 12 the full width a z in z is about half the width of 
the RF-field free case [25], whereas the width in x has 
increased by about 50%. This change of shape should in- 
fluence the ballistic expansion of the BEC once the trap- 
ping fields are switched off. Furthermore, the condensate 
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is slightly bent in the z direction. 

3. The ground state in the adiabatic potential is a 
dressed state, i.e., it is a superposition of different hyper- 
fine sublevels. Observing the different hyperfine compo- 
nents would very convincingly demonstrate the trapping 
of the BEC in the adiabatic potential. 

4. If the RF fields arc suddenly switched off the compo- 
nents Mp = 1 and 2 should perform harmonic oscillations 
in their respective magnetic trapping potentials. 

VI. TWO-DIMENSIONAL ATOM TRAPPING 
REGIME 

A. Basic considerations for 2D atom traps and 2D 

BEC 

Having introduced field-induced adiabatic potentials 
as a means to create new types of trapping potentials, 
we now turn to the question of how they can be used 
to obtain a two-dimensional atom trap by the appropri- 
ate choice of strong fields and large detunings. Unfortu- 
nately, as mentioned above, in these parameter regimes 
the resulting systems are hard to model numerically. So 
in this section we will use general arguments to show 
under what conditions it might be possible to obtain 2D 
trapping, and even a 2D BEC, using the kind of adiabatic 
loading scheme mentioned in Sec. V. These arguments 
will allow us to relax the restriction on the geometry of 
the system somewhat (the magnetic trap can now have 
three different frequencies). However, in our discussion 
we explicitly consider a two-state system (rather than 
general multi-level systems). Nevertheless we expect that 
much of the discussion of a two-state system would hold 
qualitatively for multi-level systems with an appropriate 
choice of parameters. 

In general, a 3D harmonic potential with angular 
frequencies u>i, u>2, utrans provides an effective two- 
dimensional trap for atoms of temperature T, if 

TkJi t 2 < k B T < huj trans (25) 

(from now on we return to unsealed quantities). Thus a 
good 2D trap has a large uitrans to allow one motional 
degree of freedom to be frozen out at high temperatures. 
For an ideal gas, Bose condensation in a 2D harmonic 
trap occurs at 

k B T c = hujV6N/i: (26) 

with u> — (0J1LU2) 1 / 2 and N the number of atoms [22]. 
Thus, combining Eq. (25) and Eq. (26), the number of 
atoms that can undergo a genuine 2D condensation satis- 
fies the condition TV < (w tra „,,/<I>) 2 . For a higher number 
of atoms condensation would have occurred already in 
the 3D regime. Therefore we see that, to aim for a 2D 
BEC, a high ratio w tra „ s /a) is a further desirable criterion 
for a 2D trap, especially considering that small numbers 
of atoms would be hard to image experimentally. 



B. Two-dimensional atom trapping and adiabatic 
potentials 

We will generalize the discussion in Sec. V so that we 
now allow the original magnetic trap to have different 
oscillator frequencies ui x , u> y , and u> z along the three axes 
with corresponding oscillator lengths a x , a y , a z . (Note 
that Oj = yjTi/mijJi, i = x,y,z). Thus the corresponding 
magnetic trap has the potential 

V tr (r) = ±m(ul* l +uy+ulz 2 ) (27) 

where, again, we are now working with unsealed quanti- 
ties. As in Sec. Ill, we consider a two-level system and 
assume that the second state experiences the magnetic 
potential — V tr . 

The upper field-induced adiabatic potential, which is 
the potential of greatest interest, is then given by 

V(r,t) = mgz + V+(r,t) (28) 

with 

V+ (r, t) = ^[V tr (r)-hA(t)/2}^ + (hn(t)y. (29) 

Here A denotes the unsealed detuning, i.e. A = uif — 
AE(0)/h with u) f the RF field frequency and AE(0) the 
minimum energy difference between the hyperfine states 
(see Sec. II. A). The 'bare' potentials intersect, where 
there is resonance, i.e., at the locations HA(t) — 2V4 r (r). 
This intersection we call the 'seam' of the bare poten- 
tials. In the absence of gravity, the seam of V + would 
also be the location of the minimum in V and would have 
the shape of an ellipsoid with radii r,- L = (A/cjj) 1 / 2 ^, 
i = x, y, z. Thus, without gravity, or gravity being com- 
pensated for, the atom distribution forms an ellipsoid 
bubble, or shell in the potential (29). 

Under the influence of gravity the atoms sag to the 
bottom of the shell potential (28), and this would not 
allow the formation of a matter- wave bubble. In fact, as 
shown in Sec. HID, a closed bubble could exist only up 
to radii r < (Tf/gm 2 ) 1 / 3 (~5x lCT 7 m for Rb87). Under 
typical experimental conditions we can expect a situation 
as depicted in Fig. 10 which shows — for the indicated 
parameter values — the potential V in the plane y = 
along with an atomic BEC in the ground state for this 
potential. 

Typically, an atom moving in the potential V is con- 
fined to the close proximity of the seam due to the strong 
confinement in the direction transverse to it. If its en- 
ergy energy E satisfies E -C mgr z , it is also restricted 
to the vicinity of the bottom of V because of gravity. 
To a good degree of approximation the atomic motion 
can therefore be modelled as harmonic. To determine 
the harmonic frequencies we can expand Eq. (28) about 
z = —r z to find 

Wi, 2 = {g/r z ) 1/2 u} x . y /u z (30) 
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FIG. 10: Adiabatic potential V(r) = V+(r) + mgz (in arbi- 
trary units) at y = for lo x /2-k — 11 Hz, tUy/2-K = lu z /2tx — 
220 Hz, H/27T = 2.6 kHz, and A/2tt = 66 kHz. Inserted is the 
ground state of a Rb87 condensate with about 10 5 atoms. 



for atoms held in the adiabatic trap. Equation (10) gives 
an estimate of the decay rate and if we require that the 
decay rate is less than a certain maximum, i.e., 

Tea; ^ ^fmax 

(33) 

it follows from Eq. (10), when unsealed, that one needs 

Cl 3 > Xlu 2 A (34) 

with the parameter A = 2[ln(-f max /2uj z )] 2 /n 2 . Thus the 
worst allowed decay rate can be specified in terms of the 
vertical trap frequency. In the numerical example below 
(Sec. VI D) we choose r y ma x/u z = 0.01 so that A = 5.7. 

As indicated in Sec. V.B, analytical expressions for the 
lifetime in general multi-state systems are not available. 
However, we expect that our estimates are still valid, 
at least as a first approximation, under such conditions. 
This view is supported by our numerical study of the 
five-state system. 



along the surface of the seam, and 

LO trans = (A/ft) 1 / 2 ^ (31) 

in the normal direction. Here the unsealed coupling fre- 
quency is denoted by £1. The oscillator lengths corre- 
sponding to u>i t 2 and LOtrans will be denoted by citrans, 
a\, and ai. The effect of gravity on the potential will 
also shift it downwards so that the equilibrium point is 
now at 

«>~-\F^-^— (32) 

However, the displacement of g/ui 2 rans is rather small and 
also turns out not to be sufficient to violate the harmonic 
expansion of Eq. (28) which results in Eqs. (30) and (31). 
(See section VI C 3 below.) 

C. Conditions for 2D trapping in adiabatic 
potentials 

The obvious strategy for realizing a 2D trap is to in- 
crease the detuning A as much as possible in order to ob- 
tain strong radial confinement. A sufficiently large cou- 
pling VL will be necessary to ensure that the potentials 
remain adiabatic. To arrive at a 2D trap, a number of 
considerations have to be taken into account (some given 
previously in Ref. [9]), and these factors determine how 
large Cl and A need to be, and the relationship between 
them. 



2. Strong binding and spatial thinness 

The lifetime sets a lower limit on il at a given A. 
Choosing £1 too large, however, impairs the strong radial 
confinement. An upper limit is imposed by the obvious 
conditions 

Utrans > ^1,2 (35) 

and 

T z > dtrans- (36) 

The first of these was discussed in Sec. VIA and ensures 
that one spatial direction is effectively frozen out at low 
temperatures. The second condition guarantees that the 
atoms are tightly confined to the vicinity of the seam. 
Both these requirements are typically not in conflict with 
condition (34). Formally, this can be seen from the fact 
that the first condition translates into 

- 9 h z ujl]a 2 ,„„, 
Q «^dt^> ( 37 ) 

where for u x . y one should choose the maximum of lo x and 
u>y to make the condition most restrictive. The second 
condition leads to 

Q < 4A 3 /w 2 . (38) 

Inequality (34) will allow us to choose Q oc A 1 / 3 , so that 
the above conditions are increasingly easy to meet for 
growing A, in accordance with our expectation. 



1. Lifetime 3. Harmonic 2D trapping 

For a given A, the RF coupling strength Cl has to be In this section we will look at the conditions for har- 
chosen accordingly to assure a sufficiently large lifetime monicity of the trap. Harmonicity is probably desirable, 
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but is not actually essential for trapping. The conditions 
derived below will show that the trap is operated in a 2D 
harmonic regime. 

To be able to speak of a harmonic 2D trap at a finite 
temperature, one has to require that 



hwi,2 <C k B T c -C mgr z . 



(39) 



The first criterion ensures that, at the condensation 
point, the atoms can reach sufficiently many quantum 
states to actually experience the trapping potential as 
being harmonic, ft is automatically fulfilled if we substi- 
tute for T c from Eq. (26) to find 



6 \ W 2,1 



(40) 

■Jx,y/^y,x is the ratio 



where from Eq. (30) wi, 2/^2,1 

of trapping frequencies in the x-y plane. Equation (40) 
should be easily satisfied if there are to be any significant 
number of atoms in the trap. 

The second part of Eq. (39) leads to 



5/2 

a z to z 



(41) 



when we substitute for r z . This second criterion makes 
sure that anharmonic effects in the trapping potential 
can be neglected. The same condition arises as the con- 
dition for harmonic motion of a single particle confined 
to the 'seam'. (A picture which assumes tight trans- 
verse trapping.) In the example below mgr z /kB — 18 
/zK which is much larger than T c ~ 0.3 /xK. Note, how- 
ever, that the appearance of anharmonic effects in the 
2D trapping potential does not automatically affect the 
two-dimensionality of the trap. 

We can also ask if we have harmonicity in the third, 
transverse, direction. If we focus on the bottom of the 
potential, i.e. take x = y = 0, the adiabatic potential 
(28) simplifies to 



V(z,t) = mgz + ^(mtu 2 z z 2 /2 - ftA(i)/2) + {hfl(t)) 2 . 

(42) 

Then in an adiabatic regime, one limit which produces 
quite simply a harmonic potential is found from the con- 
dition 



nn(t) > 

which together with 



1 22 

— TTLLO -Z 

2 



ftA(i)/2 



\Z - Z Q \/\Z \ < 1, 



(43) 



(44) 



allows a quadratic expansion of Eq. (42). An estimate 
for zq, the location of the minimum, has been given in 
Eq. (32). At z = z we may substitute Eq. (32) into 
Eq. (43) to eventually obtain 



Here we also assumed that the shift of z$ away from 
— r z = —\/hA/mijJ z (see Eq. (32)) was very small so 
that the smallest term could be dropped from z$. The 
condition for this approximation is that 



h~A 



muji 



> 



which would lead to the subsidary condition 

g 2 & 



A d » 



(46) 



(47) 



when we substitute for Wtrans- 

In considering values of z away from zq in Eq. (44) 
and Eq. (43), we need an estimate of the thickness of the 
matter-wave disc. A distance characterising this could be 
o-trans, i-e. we would use a single particle wave function 
width as a measure of the size of \z — zq\ in Eq. (44). 
This equation is then already satisifed by Eq. (36) and 
the remaining condition Eq. (43) is also, approximately, 
satisifed if we accept the scaling of Eq. (34), (with a 
typical value of A). 

It turns out that it is possible to satisfy all the con- 
straints in this section. It is even the case that violation 
of Eqs. (40) and (41) need not prevent a 2D BEC since 
there is no reason why a BEC cannot be 2D in an an- 
harmonic regime; equation (26) would simply not apply. 
Likewise, anharmonicity in the transverse direction, re- 
sulting in violation of conditions (45) and (47) need not 
prevent 2D trapping or a 2D BEC. 



4- Loading 

A straightforward way to load the 2D trap consists 
in starting from a condensate in the original magnetic 
potential and then adiabatically transferring it to the 2D 
trap by appropriately switching on the RF field. The 
minimum duration of this process can be estimated by 
stipulating the adiabaticity condition Co\^, <C u\ 2 . Using 
relation (34) this leads to 



t > 



Wl,2 M'y/g/OxM 



(48) 



Note that during loading, f2 can be increased above 
the intended final value to reduce intermediate adiabatic 
losses. 



5. Temperature 

For the trap to be 2D at a temperature T the second 
part of Eq. (25) applies (fc^T < Twj tra ns) which with the 
substitution of Eq. (31) leads to 



A» 



9 



(45) 



/A 



(49) 
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If we take the case of equality in Eq. (34) and substitute 
for A we find 



k B T < M/VA. 



(50) 



We would expect this condition to be more restrictive 
than Eq. (41) which leads to k B T < HQ 3 / 2 g /2VAa z w! /2 : . 
Thus, by substituting Eq. (26) into (50), the number of 
atoms that can undergo a 2D condensation is limited by 



2f)2 



6Xlu 2 ' 



(51) 



If we now utilize uj = (uji^) 1 ^ 2 and Eq. (30) we finally 
obtain 



n 2 a z 



6g\ 3 / 2 



(52) 



Since Eq. (52) puts an upper bound on the number of 
atoms it is desirable to have as large a coupling fi as 
possible. 

In fact none of the conditions above constrain f2 if the 
connection (34) is accepted. The higher the value of & the 
better (as in Eq. (52)). However, from a practical point- 
of-view, it would be desirable to find the lowest viable 
SI. This is determined by the most restrictive inequality 
on SI. That is, if we place each of our various inequali- 
ties in the form O > fc, where k contains the remaining 
constants, we will look for the inequality with the largest 
constant k. For realistic values of A and u XiV this appears 
to be Eq. (45), which on substitution of Eq. (34) leads to 



SI 7 > X 6 g 2 Lui/a 



(53) 



Any value of SI greater than this would be suitable, but 
note that increasing SI also increases the loading time, 
Eq. (48), and this is ultimately undesirable. To see what 
is realistically possible we must now determine some val- 
ues for a practical case. 



D. Numerical estimate 

We reported a numerical estimate in Ref. [9] , where we 
considered a typical Ioffe-Pritchard trap with uj x /2tt = 11 
Hz, Ldy/2ir = ui z /2ir = 220 Hz containing Rb87 atoms 
initially in the F = 2, M = 2 ground state. In order to 
match approximately the two-state theory given above 
to this five-state system we need to replace A by 4A to 
obtain the correct condition for resonance. The trap po- 
tential (27) becomes that for the Mp = 2 state, and the 
coupling SI is replaced by 2ft. Then an RF field with a fi- 
nal coupling n/2n = 15 kHz and [from (34)] A/2tt = 12.2 
MHz is one that can be provided with currently available 
technology [23]. The condition (53) is easily satisfied 
since it results in Sl/2ir ^> 405 Hz. The resulting trap 
frequencies are u> t rans/27r — 8.9 kHz, uj\/2tt = 1.3 Hz, 
W2/27T = 27 Hz. The new trap is vertically shifted by 



0.34 mm from the center of the original magnetic poten- 
tial. The critical temperature is given by 0.43/iK, so that 
up to 3.6x 10 6 atoms could be condensed. The transverse 
width of the condensate is estimated at 0.08/^m if atomic 
interactions can be neglected. The time required for the 
preparation process should be large compared to l/o»i, 
Eq. (48), i.e., of the order of several seconds. In view 
of these estimates, the experimental realization of this 
new kind of 2D trap seems to be within reach of current 
experiments. 



VII. SUMMARY AND CONCLUSIONS 

In this paper we have shown that field-induced adia- 
batic potentials can provide a robust and versatile tool to 
create novel trapping configurations for ultracold atoms. 
As specific examples, we have considered the generation 
of matter-wave shell, or bubble, states and the prepara- 
tion of two-dimensional atom traps. In a bubble state, 
the matter-wave density is localized around the surface 
of a three-dimensional sphere. We have discussed the 
preparation process which proceeds by coupling a har- 
monic magnetic trap to a repulsive harmonic potential 
by a suitably chosen time-dependent RF field. Although 
this configuration appears to be unstable at first sight, 
the bubbles are stabilized in the resulting dressed po- 
tential and their lifetime increases exponentially with 
the RF coupling strength. The bubbles can be used as 
stepping-stones for the creation of highly-excited oscil- 
lator eigenstates and 'non- linear eigenmodes.' We also 
investigated possible experiments showing collapse and 
revival effects. Although the creation of the bubbles is 
impeded by the influence of gravity and trap anisotropies, 
we have pointed out ways to overcome these difficulties 
with the help of present-day technology. 

The same principle, which has been used for making 
bubble states, can also be applied to the creation of two- 
dimensional atom traps. Under the combined influence of 
the dressing RF field and gravity, a condensate pools at 
the bottom of the resulting potential. By increasing the 
RF detuning the radial confinement becomes steeper, and 
the condensate gets increasingly squeezed until it even- 
tually reaches a 'quasi-two-dimensional' state. We have 
given estimates for the parameters necessary to reach this 
regime. 

We hope that the results presented in this article stim- 
ulate further research into the possibilities that field- 
induced adiabatic potentials offer for the creation of new 
kinds of trapping potentials and lower-dimensional ge- 
ometries. We have focussed in this paper on shell-like and 
disc-shaped traps, but by working with atomic waveg- 
uides we also expect that new tubular potentials could 
be formed. These new kinds of potentials can all be used 
to create new quantum-mechanical states of matter, and 
might also be used in the study of weakly bound clusters 
or nano-particles. 

This work was supported by the United Kingdom En- 
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APPENDIX 

In this Appendix, we outline the derivation of Eq. (8) 
for the lifetime of the bubble states. The method is de- 
scribed in detail in Ref. [14] (see also [15]) so we can 
restrict ourselves to indicating the main steps. 

Consider the time-independent linear version of Eq. (6) 



£01 



£<P2 



2cV 
1 d 2 



7T - TT <M 



A 



with e the eigenenergy. We are interested in the structure 
of the solutions at energies £ > for which resonance 
states are expected to appear. Far away from the poten- 
tial crossing at r c — \/A the eigenfunctions (f>i i2 can be 
approximated as (cf. Eq. (9) of [14]) 



ii (r < VA) ~ k\ 



b 2 (r > VA) - ft. 



Ai(oo)exp^« J k\dr + in/4 j : 
Ai(—oo)cxp^—iJ^ kidr — in/^j 

A 2 (oo)exp^i J k 2 dr — in/A 



A 2 (— oo) exp 



-/ / kidr — iir/A 
in 



+ 



(A.2) 



(2e + r 2 — A) 1 / 2 and the turning points r\ = \f& + 2e, 
r 2 = VA - 2e (or r 2 = if e > A/2). If we set, 
e.g., A 2 (— oo) equal to unity then the other coefficients 
Aj(±oo) can be regarded as scattering amplitudes for 
the interaction process inside the curve crossing region. 
Resonance states are related to poles of the scattering 
amplitudes in the complex energy plane. The determi- 
nation of the amplitudes proceeds in two steps, (i) By 
Taylor-expanding the harmonic potentials to first order 
around the crossing at r = r c , the interaction between 
the states 1 and 2 is described as a linear curve-crossing 
problem. Using the corresponding scattering matrix the 
amplitudes Ai(po) and ^(oo) can be expressed as lin- 
ear functions of A\(— 00) and A 2 (— 00). The explicit for- 
mulas (Eq. (19) of [14]) involve the quantities S and $ 
introduced in Sec. III.B. Note that the analytic form of 
the scattering matrix is only known approximately (see, 
e.g., Sec. IV of [16] for a comparison of different results), 
so that quantitative calculations based on the analytic 
approach depend somewhat on the expressions used, (ii) 
The wave function </>i has to be subjected to appropriate 
boundary conditions, i.e., <px(r = 0) = 0. This immedi- 
ately yields the relation 



Ai(oo)/Ai(-oo) = - exp[2i/3(e)] 



(A.3) 



where /3(e) = J^ 1 fci(r)dr-7r/4 = 7r(2e + A- l)/4. Com- 
bining the results of (i) and (ii) and setting A 2 {— 00) = 1 
we finally obtain (cf. Eq. (29) of [14], note the sign error) 



A 2 (oo) = 
cos /3(E) 



(A.4) 



cos /3(E) + [e 2 ^( £ ) - 1] cos^^e-^- 5 )-*^)] ' 



with the classical momenta k\ = (2e — r 2 + A) 1 / 2 , k 2 = The poles of A 2 (oo) are determined by condition (8). 
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